Spin and Charge Structure of the Surface States in Topological Insulators 
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We investigate the spin and charge densities of surface states of the three-dimensional topological 
insulator Bi^Se^, starting from the continuum description of the material [Zhang et al, Nat. Phys. 
5, 438 (2009)]. The spin structure on surfaces other than the (111) surface has additional complexity 
because of a misalignment of the contributions coming from the two sublattices of the crystal. For 
these surfaces we expect new features to be seen in the spin-resolved ARPES experiments, caused by 
a non-helical spin-polarization of electrons at the individual sublattices as well as by the interference 
of the electron waves emitted coherently from two sublattices. We also show that the position of 
the Dirac crossing in spectrum of surface states depends on the orientation of the interface. This 
leads to contact potentials and surface charge redistribution at edges between different facets of the 
crystal. 

PACS numbers: 73.20.-r, 75.70.Tj, 79.60.-i, 



I. INTRODUCTION 

A new class of materials, topological insulators, have 
attracted a great deal of attention in the last two years 
(see Refs. QUI for reviews). Typically, these are band 
insulators for which strong spin-orbit coupling leads to 
an inversion of the bulk band gap. The band inversion Q 
causes the emergence of protected two-dimensional states 
on the surface of this in other respects conventional bulk 
insulator. 

The surface states of a topological insulator have a con- 
ical energy spectrum, characteristic of eigenmodes of the 
massless Dirac equation. They are routinely described 
by the effective two-dimensional Hamiltonian 0] 





H s = v F (a x p y - CTyP X ), 



(1) 



and <j y are Pauli 



where vy is the Fermi velocity and a 
matrices. The operator <x = (a Xl a y , a z ) is usually iden- 
tified with the true electron spin [5|. The spin polariza- 
tion is then perpendicular to the electron's momentum, 
which leads to the surface states being referred to as "he- 
lical" [fj. Yet, we are not aware of a detailed investiga- 
tion of the relation between <r and the true spin. The 
nature of cr is not important for a number of phenomena 
involving the surface states, in particular to those that 
only depend on the nontrivial Berry phase acquired by 
surface-state electrons @,Hj]. However, for the properties 
that explicitly probe the electron's spin the nature of op- 
erators in the effective Hamiltonian becomes of principal 
importance 0, [To| . In this paper we elucidate how the 
spin structure of the surface electrons follows from the 
bulk Hamiltonian of the topological insulator and show 
that the relation between spin and momentum is richer 
and more delicate than suggested by a naive interpreta- 
tion of Eq. (□}. 

As a specific material we will consider the compound 
Bi2Se3, presently probably the most popular example of a 
topological insulator. This material was considered both 



FIG. 1: Polarization of a surface-state electrons on the Bi and 
Se sublattices for different orientation of the in-plane momen- 
tum p for a surface orthogonal to the quintuple layer. Spin 
textures for a symmetric model Eq. @ are shown to demon- 
strate the inevitable strong difference between the two sub- 
lattices spin for surfaces other than (111). Details of the spin 
behavior for the asymmetric case are given in Sec. IIVI 



theoretically and experimentally in Refs. [HI [HI, where it 
was shown that Bi2Se3 has a large topologically nontriv- 
ial band gap ~ 0.3 eV, leading to surface states with a 
single Dirac cone. Zhang and coworkers [12|, [l3[ also sug- 
gested a simple three-dimensional Hamiltonian describ- 
ing the long- wavelength electronic dynamics in Bi2Se3. 
This low-energy Hamiltonian contains both the true elec- 
tron spin and a pseudospin as its primary degrees of free- 
dom, where the pseudospin refers to states with support 
on the Bi and Se sublattices. 

Bi2Se3 is a strongly anisotropic material, with a lay- 
ered structure involving quintuple layers of Bi and Se 
atoms. (One quintuple layer consists of three Se layers 
strongly bonded to two Bi layers in between.) Yet, by 
virtue of their topological protection, the surface states 
exist for arbitrarily oriented crystal surfaces, not only at 
the (111) surface parallel to the quintuple layer. Based 
on the low-energy continuum Hamiltonian of Zhang et 
al. [Hj], we here show that the er-operator entering the 
surface Hamiltonian Eq. ([IJ coincides with the electron's 
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spin for the (111) surface only. For any other surface 
orientation, the spin content of the electron wave func- 
tion components is different on different sub-lattices, see 
Fig. 1. Whereas most experiments are carried out for 
the (111) surface of Bi2Se3 because the material cleaves 
well in this direction, other surfaces are also realized, for 
example in TI nanoribbons in the experiment of Ref. [7|. 
Our findings are relevant for those nonstandard surfaces. 

Experimentally, the dispersion of surface states is 
found through angular resolved photoemission spec- 
troscopy (ARPES) lH]. Spin-resolved ARPES provides 
momentum-resolved information about the spin polariza- 
tion of the surface states in topological insulators [l5l - [T7| . 
In section Hn] we consider the electron's out-of-plane po- 
larization measured via ARPES, which may be caused by 
a spin structure of the surface states that is more compli- 
cated than the spin structure suggested by Eq. (Q}. The 
out-of-plane polarization may arise from the interference 
effects of photoelectrons emitted from the two sublat- 
tices for surfaces other than the standard (111) surface, 
or from photoelectrons emitted from Se atoms for sur- 
faces that are neither parallel nor perpendicular to the 
plane of the quintuple layers." 

Another feature of realistic three-dimensional topolog- 
ical insulators, that is easily overlooked in the effective 
surface Hamiltonian ([1} (but present in the effective low- 
energy Hamiltonian of Ref. Il2l ) , is that the combination 
of broken particle-hole symmetry and angular anisotropy 
of the bulk Hamiltonian leads to different energies of the 
Dirac points at different facets of the topological insula- 
tor crystal. Different positions of the Dirac crossing (or 
the neutrality point) with respect to Fermi energy would 
result in different surface electrons densities at differ- 
ent surfaces of the crystal. However, an uncompensated 
charge density at one or more crystal surfaces comes at 
a large electrostatic energetic cost, and one expects a 
charge redistribution in order to reduce the Coulomb in- 
teraction energy. Below we describe the solution of the 
corresponding electrostatic problem at an edge of a crys- 
tal, at the intersection of two differently oriented surfaces. 
For short facets, where the Coulomb interaction energies 
are not large enough to induce full screening, a finite car- 
rier density may remain, illustrating the impossibility to 
tune the full topological insulator surface to the Dirac 
point by means of doping or gating. 

II. SURFACE STATES 

We first discuss a simplified version of the effective low- 
energy continuum Hamiltonian of Ref. fl2l which consists 
of a three-dimensional Dirac-likc Hamiltonian with mo- 
mentum dependent mass, 

H = t z M(K) +t x ((t ■ K)A, M(K) = M - BK 2 , (2) 

where we use the capital K = (K x , K y , K z ) for the three- 
dimensional momentum. The z axis is chosen in the (111) 
direction, i.e., perpendicular to the plane of the quintuple 



layers. The vectors er = (er x , a y , <x z ) and r = (t x ,t v ,t z ) 
contain two sets of Pauli matrices operating in different 
pseudospin spaces. The matrices 7V, Ty, and t z refer to 
the P1+ and P2~ orbitals of Refs.[ll[ll which are such 
that states with t z = 1 have support mostly on the Bi 
sublattice and states with r 2 = — 1 have support mostly 
on the Se sublattice. Below we will refer to these t z = 1 
and t z = — 1 states simply as states on the Bi and Se 
sublattices, respectively. The electron's spin is expressed 
in terms of the Pauli matrices as [r| 

°z — 2 ' x ~ 2 ' _ 2 ' 

(The correct form of the spin operator for Bi2Se3 is miss- 
ing in Ref. Hit ) Equation (O implies that the operators 
er/2 and s coincide for the electron residing on the Bi 
sublattice, whereas the relation between s and cr/2 in- 
volves a rotation by a it angle around the z axis on the 
Se sublattice. 

The Hamiltonian of Eq. ^ has a higher symmetry 
than the true low-energy Hamiltonian of Bi2Se3: It is not 
only valence-band-conduction-band (particle-hole) sym- 
metric, but also symmetric under spatial rotations, pro- 
vided the Pauli matrices er are rotated together with the 
coordinates, whereas the Pauli matrices r are kept fixed. 
The rotation symmetry is a mathematical artifact of the 
pseudospin basis used in Eq. ©. The strong asymmetry 
of the real crystal of Bi2Se3, in which the z axis plays a 
special role, is reflected in the anisotropic assignment of 
the spin operators s x , s y , and s z in Eq. ([3]). Corrections 
to this simplified model Hamiltonian will be discussed in 
Sec. IIVI Note that, a direct association of the spin with 
er in Eq. ([2]) would violate the parity of the true system. 
That is why one needs to introduce the asymmetric spin- 
operator with the special axis z, Eq. even in case of 
the most symmetric model Hamiltonian [l3|, [l8[ . 

The Hamiltonian Eq. ([2]) is of the second order in spa- 
tial derivatives (momentum). Following Ref. [l2l. as well 
as a number of subsequent articles [j| [l(| Il9l. |20|| we re- 
quire that all four components of the wavefunction van- 
ish at the crystal surface. This choice of the boundary 
conditions guarantees the existence of a branch of Dirac- 
like surface states with Dirac crossing at the T-point 
(the point of vanishing in-plane momentum) , even in the 
strongly asymmetric case of Sec. IIVI We should point out, 
however, that this boundary condition is not unique, and 
that other choices of boundary conditions have been ad- 
vocated in the literature. References [l8l l2ll. I22I suggest to 
use an effective Hamiltonian that is linear in momentum, 
with boundary conditions for which only one pseudospin 
component of the wave-function vanishes at the surface. 
The two choices of boundary conditions agree qualita- 
tively for the standard termination at the (111) surface 
(both give a Dirac-like branch of surface states near the 
r point), but it is not obvious how to extend the bound- 
ary condition of Refs. 0, [2l|, [22| to surfaces of arbitrary 
orientation, and the two boundary conditions may give 
different predictions for other surfaces or non-standard 
surface terminations |22j. 
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For the explicit calculation of the surface states from 
the model Hamiltonian @ it is convenient to rotate the 
coordinate system, such that the TI fills the half-space 
z < 0. This rotation does not change Eq. @, but it 
changes the relation between the spin operators and the 
cr matrices: In the rotated coordinate system, Eq. ([3]) 
takes the form 



(cr • n)n a — {a • n)n 



(4) 



where n is the unit vector pointing in the direction per- 
pendicular to the quintuple layer plane in the rotated 
coordinate frame. 

Let us introduce two two-component spinor functions 
ipB and V'S that describe the pseudospin r z = 1 and 
t 2 = — 1 components of the 4-component wave function. 
Eigenfunctions of the Hamiltonian Eq. ^ at energy e 
that correspond to surface states with momentum p par- 
allel to the surface should be found as a linear combina- 
tions of the functions 

4) B = u B e i( - p * x+p y y}+Xz , ip s = u s e i{p * x+p v v)+Xz , (5) 

with different values of A (Re A > 0). Here and below, 
we use the lower-case symbol p = (p x ,P y ,0) to denote 
the in-plane momentum of the surface states. The spinor 
amplitudes, ub and us, satisfy the coupled system of 
equations 

(e-M(p,X))u B = A(a x p x + <7yPy - iXa z )u s , 

(e + M{p, X))us = A(a x p x + a y p y -i\a z )u B , (6) 

where M{p, A) = M - Bp 2 + BX 2 and p 2 =p 2 x +p 2 y . The 
inverse decay length can take two values, A1.2, which are 
the solutions of the equation 

e 2 = (M - Bp 2 + BX 2 ) 2 + A 2 p 2 - A 2 \ 2 . (7) 

with Re A > 0. In addition, the surface state satisfies 
the boundary condition iPb{z = 0) = ips{ z = 0) = 0, 
which provides one additional constraint from which the 
dispersion e(p) can be calculated. 

To find the solution of Eqs. (JOJ) and (O one may first 
guess the (correct) result for the spectrum, 

e(p) = ±A P . (8) 

Then Eq. ([7]) yields a quadratic equation for A, 

AX = M( P , A) = (M - Bp 2 + BX 2 ), (9) 

where the parameters A and B are taken from the orig- 
inal Hamiltonian Eq. ©. The first equality here is the 
energy of the surface state consistent with the surface 
Hamiltonian Eq. (fTJ) with A — vp. The second equality 
gives the values of 



Ai, 2 = (A± V 'A 2 + 4B 2 p 2 - AMB)/2B, 



(10) 



consistent with this value of energy. Substitution of the 
ansatz ([5]) into Eq. ([5]) then gives identical spinors ub and 



us for the two roots A1.2, which guarantees that an ap- 
propriate linear combination can be found that satisfies 
the boundary condition 4>b{z — 0) = ips(z — 0) = for 
all four components of the spinor wavefunctions ips and 
ips simultaneously. The spinor structure corresponding 
to the cases e — ±Ap is 



1 



u B = 



±i 
7 



us = 



=R 

7 



(11) 



with 7 = (p x 
the form 



ip y )/p. The full surface state then has 



ipB,s(*) = u B ,se 



)• (12) 



For a particle-hole symmetric Hamiltonian Eq. ^ an 
electron in a surface state is found with equal probability 
on either Bi or Se sublattice. In order to clarify the spin- 
content of the surface states, we first calculate the expec- 
tation value of the in-plane components of er-operator on 
each sublattice. For the in-plane components one finds 



(<Ti) B = U B (JiU B = ± 



zPj 



(a t ) s = u' s aiUs 



2p 

Vi)B, 



where i. j — 



J = x,y, 



whereas 

(<Jz)b = {&z}i 



0. 



(13) 



(14) 



As in Eq. (fTT|) . the upper/lower sign corresponds to the 
energies above/below the Dirac crossing. We see that 
for both sublattices the expectation value of er is an in- 
surface-plane vector perpendicular to the momentum p. 
However the direction of (cr) is opposite for two sublat- 
tices. As we see from Eqs. © and ((!]), the direction 
of (<t)b for the Bi sublattice always coincides with the 
direction of the true spin. Thus the spin of this com- 
ponent is always described by the surface Hamiltonian 
Eq. (fTJ. On the contrary, the spin of Se sublattice elec- 
tron component of the wave function differs from ((Ji)s 
by a 180-degrees rotation around n, the axis normal to 
the layer plane. 

We conclude that the spin orientation of two sublat- 
tices coincide with each other and with the prediction of 
Eq. ([l]) if and only if the surface of the crystal coincides 
with the quintuple layer plane. For any other crystal sur- 
face the two sublattices spins differ, and the difference is 
typically large. In Fig. 1 we show the spin-directions for 
two sublattices for a surface normal to the layer plane. 

For certain applications, such as the description of spin 
and angle- resolved photoemission (see next Section), it is 
preferable to write the Hamiltonian ^ in a form in which 
the electron spin operator s is directly proportional to 
cr. (Photo-electrons do not carry a sublattice index, so 
that a spin-operator that contains r is problematic in 
that context.) Hereto, following Ref. [1J|, we perform 
the unitary transformation 



H = UHU*, U = ^£ + i(a . n)i-^ 



(15) 
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where n is the vector normal to the quintuple layer. The 
Hamiltonian H no longer has the manifest rotational 
symmetry of Eq. ([2j , but the relation between spin and 
Pauli matrices now takes the standard form 



s = <x/2. 



(16) 



After the unitary transformation the two spinor ampli- 
tydes take the form 

u B = u B , u s = i(a ■ n)u s , (17) 

where ub and us are given in Eq. (Ilip . 

III. PHOTO-EMISSION 

Spin and angle resolved photoemission spectroscopy 
(spin-resolved ARPES) [23| has been the tool of choice 
for the experimental investigation of the helical nature of 
the surface states in three-dimensional topological insula- 
tors [H, HH . What does the coexistence of two different 
spin polarizations of the surface electron imply for the 
polarization of photoelectrons? 

In general, the probability P of electron photoemis- 
sion is proportional to the squared matrix element of the 
interaction H mt with the photon field, 



P oc 



(18) 



where |f) and |i) are the spinor states corresponding to 
the free final and the initial surface electron states, re- 
spectively. (One has H mt — (e/c)A • j, where A is the 
vector potential of the photon field and j is the nonrel- 
ativistic current density operator.) In our case, the ini- 
tial state is described by separate amplitudes ub and us 
for electrons from the Bi and Se sublattices, which are 
given in Eqs. (|11I17[) above. Since the interaction H mt is 
spin-conserving and local, the spinor structure Uf of the 
photoelectron becomes a linear combination 



U{ — au B + pus 



(19) 



where the coefficients a and ft contain contributions from 
matrix elements of H int on the Bi and Se atoms in the 
crystal, respectively. Their precise value depends on the 
details of the atomic structure and the frequency of the 
incident light. 

Since the surface Hamiltonian Eq. ([1} predicts a spin 
polarization parallel to the surface, direct observation of 
an out-of-plane spin would be of most interest. A simple 
calculation gives for the z-component 



ufa z Uf 



(n z \ft | 2 + Reaft*)(n x p y - n y p x ) 



P 



(20) 



We see that in case of arbitrary surface orientation the 
electron emitted from the Se sublattice acquires an out- 
of-plane polarization even without interference. The na- 
ture of this out-of-plane polarization is a direct conse- 
quence of the special form of the spin operators in Eq. Q : 



The surface electrons described by the symmetric Hamil- 
tonian Eq. for each sublattice have the expectation 
value of the er-operator parallel to the surface. However, 
since the spin of one sublattice is obtained from <x/2 by 
a 7r rotation around the axis normal to the quintuple 
layer, n, the spins of this sublattice all lie in a plane dif- 
ferent from the crystal surface plane. (Similarly for the 
asymmetric Hamiltonian of the following section, for the 
general surface orientation, spins of both sublattices at 
different values of momentum lie in two planes, different 
from each other and from the surface plane.) 

As we saw in Sec. HU only for the most studied (111) 
surface (the surface parallel to the quintuple layer), po- 
larizations of the electron on both sublattices coincide. 
Consequently Eq. (|20[) predicts no out-of-plane polariza- 
tion here. 

Most interesting is the photo-emission from the sur- 
face normal to the quintuple layers (n z = 0), where 
the out-of-plane polarization appears due to the inter- 
ference of contributions from the two sublattices only, 
ufa z Ui ~ Re aft* (j2"Uj) . The unit cell of Bi 2 Se3 consists 
of 5 atoms, each hosting a p z electron. The Hamiltonian 
Eq. operates in a reduced 2x2 pseudospin subspace 
built from an even (with respect to reflection through 
the middle plane of the quintuple layer) state from the 
conduction band and an odd state from the valence band 
12]. In our case of a surface normal to the quintuple 
layer, in the odd state electron waves emitted from the 
different atoms would cancel each other, unless there is a 
finite phase difference due to the momentum of outgoing 
electron. Thus we expect the photoionization amplitude 
from the odd state to be ft ~ (p • n). Consequently the 
out-of-plane spin component Eq. (f2T3f should have a node 
and change sign at (p • n) =0. 



IV. ANISOTROPIC TOPOLOGICAL 
INSULATORS 

In the simplified model of Sec. [Hi the origin of the 
different spin content of the surface states at differ- 
ent facets was the special form of the spin operator 
Eq. (U]). Whereas the special symmetries of the Hamilto- 
nian Eq. ([2]), rotational and particle- hole symmetry, are 
broken in real crystal, the mechanism leading to the dif- 
ferent spin content of surface states at different facets 
continues to operate, as we now show. 

The general anisotropic and particle-hole asymmetric 
generalization of the Hamiltonian © is [24 1 



H = e(K)+T z M+T, x J2A l a l K l , (21) 



x,y,z 



where now 



and 



M(K) =M-^2 KjSjjKj 



(22) 



(23) 
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with a positive-definite symmetric matrix Bij and a sym- where 



metric matrix D^. The first term e(K) in Eq. ([2~Tj) is 
explicitly particle-hole (valence band-conduction band) 
asymmetric. Choosing the z axis to be orthogonal to the 
plane of the quintuple layer of the Bi2Se3 crystal, the 
matrices Bij and Dij are found to be diagonal fl2l ] 



Bij Bi&ij , Dij DiSij : 



(24) 



The microscopic calculations of Ref. [12| give the follow- 
ing values for their elements: 



(25) 



A z 


= Ai 


= 2.2eVA, 


B z 


= Bx 


= lO.eVA 2 , 


D z 


= D 1 


= 1.3eVA 2 



for the direction normal to the layer, 

A x = A y = A 2 = 4.1eVA, 
B v = B 2 = 56.6eVA 2 



B.i 
D, 



(26) 



D y = D 2 = 19.6eVA 



for the two axes in the layer plane, and 
AI = 0.28eV. 



(27) 



Finding the surface states for the general Hamiltonian 
Eq. (|2"Tj) become a complicated algebraic problem. Here 
we discuss the solutions only in case that and D^ 
remain diagonal after a rotation that brings the surface 
normal to the z axis. This includes the case of a surface 
parallel to the quintuple layer and a surface perpendic- 
ular to the layers in E^Sc;), as in the nanoribbons of 
Ref. 0. 

As in SecHH we search for eigenfunctions of the Hamil- 
tonian (|2"Tj) of the form (we use again K for the three- 
dimensional momentum and p for the two-dimensional 
surface states momentum) 



ti'i 



u B e 



i(p x x+p y y)+\z 



use 



i{p :L x+p y y) + \ Z 



(28) 



where Re A > and we rotated the coordinate system, 
such that the topological insulator occupies the half space 
z < 0. With the particle- hole asymmetric term e(K), 
simple guessing of the appropriate solution, like we did 
in Sec. ITU does not work and one has to pursue an explicit 
derivation of the result. To this end it is convenient to 
replace the sublattice spinors us — (MB t ,«Bj and us = 
(ws T , us l ) by spin- up and spin-down pseudospinors 



u = (u Bf ,u Sf ), v = (u Bl ,u Sl )- 



(29) 



These pseudospinor amplitudes now satisfy the system of 
equations 

(e - e - Mt z + i\A z t x )u = (A x p x - iA y p y )T x v, (30) 
(s-e- Mt z - i\A z r x )v = (A x p x + iA y p y )T x u, (31) 



M = M(p,X)=M -B xPx \~Byp 2 y + B z \ 2 , (32) 



From the first equation we find 

v = Qu 

with 

1 



D x pl + Dyp 2 y - D Z X 2 . 



n = 



(e - e)r x + iMTy — i\A z 



(33) 



(34) 



A X p X iAyPy 

Substituting this expression for v into Eq. pip one finds 

[e - e(p, A)] 2 - M( P , A) 2 + A 2 p 2 + A 2 yP 2 y - A 2 \ 2 . (35) 

This equation is biquadratic in A and allows us to find 
two values of the squared inverse decay length A 2 2 for 
each energy e, and consequently two decaying solutions 
with Re Ai,2 > 0. 

The open boundary conditions at the surface can be 
satisfied if there exists a choice of pseudospinors corre- 
sponding to two solutions Ai and A2 such that 



Ui = u 2 , Vi 



v 2 - 



(36) 



The first equality here is easily satisfied, because Eq. ((35]) 
imposes no restrictions on the pseudospinor u. The sec- 
ond equality then involves the fi = f2(e, A) of Eq. (|34p . 
which depends on energy and on the specific root of the 
biquadratic equation (|35p . 

(37) 



= ui-«2 = [n(e,Ai)-fi(E,Aa) 



This equation has a solution only if the difference ma- 
trix fi(e, Ai) — fi(e, A2) has a zero eigenvalue, and con- 
sequently det[f2(e, Ai) — fl(e, A2)] = 0. Thus we find the 
condition 

det[(e 1 -e 2 )T x + i(Mi-M 2 )T2-iA z {\ 1 -\ 2 )] =0, (38) 



where Ai and A2 are the two roots of Eq. (f3"5]l . This 
immediately gives a simple formula 



A 2 z = {B 2 z -D 2 z ){\ 1+ \ 2 f 



(39) 



After straightforward calculation we can now find the 
energy of the surface states 



e(p) 



D 



Dl 



B M±^A 2 x p 2 x + A 2 yP 2 y Jl - B2 



B z 



D Z B X 2 



D y B z 



B 



(40) 



and the two spinor amplitudes (here we use again the 
more meaningful spinor functions for two sublattices in- 
stead of the pseudospinors u and v) 



UB = 2V 1+ 57 



1 A D * 
us = 2T-B: 



±1 

=R 

1A 



(41) 
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where now 



1A 



^ {A X P X + iAyPy)/(A X P X - lAyPy). (42) 



The main differences between the solution of the sim- 
plified Hamiltonian of Sec. [TT] and the solution given in 
Eqs. (gnj) and gj) are: 

(i) The electron no longer spends equal time on the two 
sublattices. The difference in probabilities is governed by 
the ratio D z /B z , which can actually be small. (Reference 
GJ finds Di/Bx w 0.13.) 

(ii) For the anisotropic Hamiltonian the expectation 
value of the er-operator, which we discussed in section [TT1 
is no longer perpendicular to the momentum. However, 
the expectation values {<t)b,s on the Bi and Se sublat- 
tices are still perpendicular to the rescaled momentum 
p x — > A x p x , p y — > AyPy, and Fig. Q] remains valid if 
the rescaled momentum is used. We observe that such a 
rescaling can also be used to account for the difference 
between the surface-state spin on the Bi and Se sublat- 
tices that is found here and in Sec. [II] Rescaling with one 
inverted sign, p x — > —A x p x , p y — > A y p y , is sufficient to 
permute the two panels in Fig. [TJ 

(Hi) The third important feature of the solution 
Eq. (|4"0|) is also caused by the particle-hole asymmetric 
term eo(K) in Eq. (|2"Tj) . As long as one considers the 
bulk states, this term leads to a trivial bending of both 
conduction and valence bands. However, the inclusion 
of the energy eo(K) into the Hamiltonian (f2~Tj) leads to a 
nonzero value of the Dirac crossing energy, as can be seen 
from the explicit solution (|40|. In the most general case 
of matrices Dij and Bij having no vanishing elements 
we found the position of the Dirac crossing, which is the 
(unique) energy eigenvalue of a surface state with van- 
ishing momentum parallel to the surface, to be simply 
e — EDirac = MD ZZ /B ZZ . Thus, returning to the original 
coordinate system of Eqs. ([2"Tj) - ([2"T|) , we may introduce a 
vector e normal to the surface at a particular point and 
write a formula for the energy of the Dirac crossing valid 
for the entire crystal 

EDirac = M C^2 e i D ij e o) I '(^ <, . / '. . <> ; : - (43) 

Substituting the actual values of the matrices and 
Bij, we find that for Bi2Se3 the position of the Dirac 
crossing may vary by as much as 



^_^V«.06eV«^-, (44) 
B 2 Bij 6.8nm' v ' 



upon varying the orientation of the surface. 



V. ELECTROSTATICS OF A TOPOLOGICAL 
INSULATOR EDGE 

Being built of neutral atoms, the topological insulator 
crystal is obviously a charge-neutral system. Because of 
the bulk band gap, the bulk charge density is insensitive 




FIG. 2: Electric field near the edge of an anisotropic TI. 



to the precise choice of the Fermi energy. However, the 
requirement of neutrality should also include the surface 
charges. Even in the absence of particle-hole symmetry, 
the structure of the surface-state spectrum is such, that 
generically charge neutrality is achieved if the Fermi en- 
ergy equals the energy of the Dirac crossing, EDirac- Since 
the latter energy depends on the orientation of the sur- 
face, see Eqs. ([4"3l and ([SI) , surface-charge neutrality 
can not be satisfied by a uniform choice of the Fermi en- 
ergy for the entire crystal. (The variation of the surface 
charge density for a hypothetical uniform choice of the 
Fermi energy is estimated as An ~ (AeDirac/ A 2 ) 2 / Air w 
1.7 • 10 n cm~ 2 .) Such a choice of the Fermi energy would 
lead to a strong electric field along the surface, which 
should not happen for a metallic surface. Instead, the 
surface charge is redistributed near the edges between 
different facets, so that the electric field parallel to the 
metallic surface vanishes. Far away from the edge, sur- 
face neutrality will then be achieved simultaneously for 
both surfaces, while the electrostatic potential difference 
between two facets will compensate the difference of the 
Dirac crossing energies, A£rji rac = AV. In the limit that 
the surfaces can be approximated as perfect metals, the 
solution of such electrostatic problem is straightforward 
25]. For two surfaces at a 90 degree angle, as shown 
in Fig. O the potential 3>i n and $ ou t inside (polar angle 
— 7r/2 < <fr < 0) or outside (polar angle < 4> < 3n/2) 
the topological insulator is 



$out(r, 



2AV 

Sire 



0, $in(r, 



2AV 



corresponding to a surface charge density 
i / m / I, AV 



3 2~K 2 er 



(45) 



(46) 



at a distance r away from the edge. Substituting the 
dielectric constant of Bi2Se3, k w 110, and taking 
AV = .06eV from Eq. (|44[) . we arrive at the estimate 
\<j(r)\ ~ (e/r) x 2.3 • 10 6 cm _1 . Doping of the surface elec- 
tron gas by electrostatic external gates would lead to an 
additional edge/corner charge accumulation, similar to 
the edge charge accumulation in graphene ribbons (26[. 

The electrostatic calculation of Eqs. (I45I46P ignores 
completely the kinetic energy of the surface electrons. 
However, due to the large dielectric constant in Bi2Se3 
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there exists a regime where the semiclassical electron dy- 
namics is already described by a smooth coordinate de- 
pendent Fermi energy Ep(r) — Akp(r), but the electro- 
static energy is not yet dominant, e$ ~ Ep (Thomas 
Fermi approximation). In this case one requires the 
sum Ef{t) + e$ to be a constant over the metallic sur- 
face, while electrostatics is recovered in the large-sample 
limit. The range of validity of Eq. (1451) is now found as 
AV > A 2 k F (r) = A 2 2y/na(r)/e, leading to r > 135nm. 
For much smaller samples one should ignore the Coulomb 
interaction and describe the surface electrons by the 
Dirac equation with shifted crossing energy, Eqs. (|43|44[) . 



VI. DISCUSSION 

In this article, we have analyzed the spin structure of 
surface states in the three-dimensional topological insu- 
lator Bi2Se3, as it follows from the low-energy continuum 
description proposed in Ref. [H The inequivalence of the 
spin variables on the Bi and Se sublattices in this con- 
tinuum description leads to a nontrivial spin structure of 
the surface states at surfaces other than the (111) surface. 
We also found that the energy of the Dirac crossing in 
the surface-state dispersion depends on the surface ori- 
entation. Although the precise boundary conditions of 
the effective low-energy description — open boundary 
conditions, in which all components of the spinor wave- 
function vanish at the surface of the topological insula- 
tor, as in Ref. — are key to our quantitative analy- 
sis, we expect that the effects we predict persists if the 
boundary conditions are changed. Different boundary 
conditions hav e ap peared in the literature for the (111) 
surface but not for other surfaces of a Bi2Se3 

crystal. 

The theoretical findings of this paper may be applied 
to the interpretation of several experiments involving the 
surface states of three-dimensional topological insulators. 
Here we mention the STM measurements of the surface 
electrons charge near an artificial step on Bi2Te3 sur- 
face [13, [H| and the measurement of Aharonov-Bohm 
oscillations in the surface-state-mediated transport in 
Bi2Se3 nanoribbons @- In the latter case a theoretical 
analysis predicts that the magnetoconductance should 
depend strongly on the position of the Dirac crossing for 
the electrons on the surface of the nanoribbon Q . Obvi- 
ously, a non-uniformity of this Dirac crossing of the type 
discussed above will affect the minimal conductivity and 
interference pattern for the thick (~ lOOnm) rectangular 
shaped nanoribbons of Ref. Q and should be taken into 



account in a quantitative modeling of the device. 

Our most remarkable result is the prediction of the pos- 
sibility of out-of-plane momentum polarizations of photo- 
emitted electrons, Eq. (|20|) . depending on the exact orien- 
tation of the surface. Measuring these spin components 
would confirm the validity of the microscopic Hamilto- 
nian of Ref. 12[ and its boundary conditions. It will 



also explicitly demonstrate the signatures of interference 
of photoemission contributions from two sublattices of 
Bi2Se3. Such experiments require, however, the prepa- 
ration of Bi2Se3 crystals with sufficient quality surfaces 
other than in the (111) direction. 

Although the low-energy model of Ref. El that we used 
as the basis for our calculations, predicts a nontrivial spin 
structure for surfaces other than the (111) surface, there 
may be other mechanisms that lead to an effectively re- 
duced spin that are not included in this model. In this 
context, we mention recent first-principles calculations of 
Yazyev and coworkers, who find that the net spin polar- 
ization of surface states on the (111) surface is reduced 
by an amount of order 50% [2!| , which they attribute to 
the effect of strong spin-orbit interaction on the heavy 
Bi atoms. Since the magnitude of the spin of a photo- 
emitted electron has to be s = 1/2 (a free electron always 
has its spin pointed in some direction), the reduced net 
spin found in Ref. [2^ implies a likely out-of-plane compo- 
nent of photo-emitted electrons, even for the (111) sur- 
face. This effect may be suppressed since it involves the 
photo-ionization from the inner row of Bi atoms, less ac- 
cessible for ARPES. On the other hand, as the spin-orbit 
interaction admixes p x ± ip y electron states to p z ones, 
one may play with the photon polarization to selectively 
enhance the electron's ionization from Bi. Further inves- 
tigation in this direction is obviously desired. 
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